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ABSTRACT 



Rate coefficients for inelastic Na+H collisions are calculated for all transitions between the ten levels up to and including the ionic 
state (ion-pair production), namely Na(3s,3p,4s,3d,4p,5s,4d,4f,5p)+H(ls) and Na + +H~. The calculations are based on recent full 
quantum scattering cross-section calculations. The data are needed for non-LTE applications in cool astrophysical environments, 
especially cool stellar atmospheres, and are presented for a temperature range of 500-8000 K. From consideration of the sensitivity 
of the cross-sections to input quantum chemical data and the results of different methods for the scattering calculations, a measure of 
the possible uncertainties in the rate coefficients is estimated. 
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1. Introduction 

Lack of accurate data for collision processes needed for re- 
liable non-LTE line formation calculations in cool star atmo- 
spheres, especially processes involving hydrogen atom and elec- 
tron impacts, poses a m ajor sourc e of un c ertainty fo r stella r 
abundance ana lyses; e.g. Lambert (1993); iKiselmanl (1200 ll) ; 
lAsplundl (12005). For hydrogen, the situation is particularly poor. 
The possible importance of collisions with neutral hydrogen in 
non-LTE line formation calcu lations was first pointed out by 
ISteenbock & Holwegerl (1 19841) in their study of the of the statis- 
tical equilibrium of Li in cool stars. Although inelastic processes 
due to hydrogen collisions are expected to be much less efficient 
than those due to electrons, the much greater abundance of hy- 
drogen atoms may overcome this reduced efficiency: in the line 
forming regions of a solar-type star the abundance of hydrogen 
atoms is typically four orders of magnitude greater than that of 
electrons, even more in metal-poor stars. 

At the time of Steenbock and Holweger's study there was 
practically no reliable experimental or theoretical work on in- 
elastic hydrogen collision processes; however, the situation has 
been improving slowly but steadily over the quarter of a cen- 
tury since then. T heir work prompted an experimental study by 
iFleck et al.l (1199 lh of Na(3s) + H -> Na(3 p) + H at low ( 1 5-1 500 
eV) energies, though due to experimental difficulties not down 
to near the threshold (2.1 eV for this case), which is the relevant 
regime for the temperatures of interest in cool stars. Revised ex- 
p erimental data, i nclud ing results down to 10 eV, were presented 
in lBelvaev et al.l (1 19991) . This work was followed by a number of 
theoretical studies involving some of the present authors. First, 
quantum scattering calculations were performed for Na(3s) + 
H - > Na(3/7,4s) + H down to the threshold (Belyaev et al. 
1999) and found good agreement with the experimental results 



above 10 eV. T his work on Na+H was f ollowed by calcula- 
tions for Li+H dBelvaev & Barklemll2003[) based on qua ntum- 
chemical data calculated by some of us (Croft et al. 1999) . This 
was followed b y astrophysical application (Barklem et al, 2003; 
Lind et al . 2009), where it was found that excitation collisions 
Li(nZ) + H — > Li(«'/') + H were unimportant, yet the ion-pair pro- 
duction and mutual-neutralisation process Li(3s)+H ^ Li + +H 
(often referred to as charge exchange in astrophysics) was found 
to be rather important, resulting in changes in spectral line 
strengths of around 20% in cool, metal-poor , sub-giant stars. 

In a recent paper (Bel vaev et al.l l20Toh we revisited low- 
energy Na+H collisions, since for astrophysical non-LTE mod- 
elling of Na line formation data for transitions between all pos- 
sible Na levels are needed, while the earlier experimental and 
theoretical studies dealt primarily with the resonance transition 
(which corresponds to the Na D lines). Details of the calcula- 
tions can be found in that paper. Cross-sections for transitions 
between all ten levels up to and including the ionic state (ion-pair 
production) for collision energies from threshold to 10 eV were 
presented. In fact, cross-sections were calculated up to collision 
energies of 100 eV; however, the results at energies higher than 
10 eV are of little importance at the temperatures of interest. The 
purpose of this research note is to present rate coefficients cal- 
culated from these cross-sections, as these rate coefficients are 
needed for non-LTE applications such as in cool stars. 



2. Collision rate coefficients 

The rate coefficients, (crv), for excitation and deexcitation pro- 
cesses Na(«Z) + H(ls) ^ Na(n'Z') + H(ls), and for the ion- 
pair production and mutual-neutralisation processes involving 
the ionic state Na(n/) + H(ls) ^ Na + (2s 2 2/? 6 ) + H , are pre- 
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sented in TableQ] The coefficients have been obtained by folding 
the cross-sections with a Maxwellian velocity distribution from 
threshold to 100 eV, and are presented for temperatures in the 
range 500-8000 K. 

The 500 K data are provided since there is substantial interest 
in Na lines in low temperature astrophysical envi ronments such 
as brown dwarfs and plan etary atmospheres (e.g. lBurrows et all 
1200 it iBarman et al.ll2002h . and though other perturbers such as 
H2 and He are usually more ab undant in these situations, the data 
may be useful. lLin et alJ d2008l) have recently calculated data for 
inelastic processes in Na+He collisions with application to plan- 
etary and brown dwarf atmospheres in mind. However, the main 
driver behind our study is the need for data for line formation 
modelling in F, G and K star atmospheres where ground state hy- 
drogen atoms are indeed the most abundant perturber. Thus, our 
discussion below will predominantly focus on the 2000-8000 K 
data. 

Although the data for inverse processes are related via the 
detailed-balance relation and thus redundant, little space can 
be saved by omitting one of the two and so we present data 
for both, thus saving the end user the need to make such cal- 
culations themselves. It is worth noting that the excitation and 
ion-pair production processes were calculated directly from the 
cross-sections, while dexcitation and mutual-neutralisation pro- 
cesses were calculated via the detailed-balance relation rather 
than di rect calculation from the deexcitation cross-sections. As 
for Li dBarklem et al.ll2003h . for the temperatures most relevant 
for cool stars around 6000 K, it is readily seen that the rate coef- 
ficients for ion-pair production are remarkably large, particularly 
from the first excited s state, here 4s. Some rates between close- 
together neighbouring levels, such as 4s — > 3d and Ad — > 4/ are 
also seen to be very large. At low temperatures, the rate coeffi- 
cients for these transitions with small energy thresholds become 
by far the largest. 

In iBelvaev et alJ (l2010h it was shown that narrow orbiting 
and Feshbach resonances are often present in the cross-sections 
at low energy (see Figs. 5, 8 and 9 of that paper). To fully resolve 
these features requires calculations with very small steps in col- 
lision energy, which is extremely time consuming. Such detailed 
calculations were performed for the 3 s — > 3p cross-sections in 
order to resolve a significant fraction of the orbiting resonances. 
This further allowed pure background cross-sections, i.e. with no 
resonances, to be constructed for this transition. From these cal- 
culations it was estimated that the contribution of the resonances 
to the rate coefficient was roughly 33% at 500 K, 20% at 2000 K, 
15% at 5000 K and 10% at 8000 K. The rate coefficients for 
3 s — » 3p presented here are calculated including the contribution 
of orbiting resonances. In all other cases, the cross-section cal- 
culations were performed using an energy grid with steps much 
greater than the widths of the resonances. Consequently most 
resonances are missed and those resonances that are by chance 
sampled are not resolved. Missed resonances will lead to an un- 
derestimate of the contribution of the rate coefficient, while un- 
resolved resoances will generally lead to an overestimate of the 
contribution since the widths of the resonances are significantly 
overestimated, and this may lead to some cancellation of errors. 
In any case, the error in the rate coefficients from this source of 
uncertainty should not exceed the magnitude of the corrections 
given above for 3s — > 3p, i.e. of order 10%, which is signifi- 
cantly less than the uncertainties from other sources discussed 
below. 

The accuracy of these rate coefficients is determined by the 
accuracy of the cross-sections near the threshold. At T = 2000 K 
the values of the rate coefficients vary from 10 -31 to 10~ 9 cm 3 /s, 



more than 20 orders of magnitude. At T = 8000 K the varia- 
tion is from 10~ 22 to 10~ 9 cm 3 /s, 13 orders of magnitude. This 
reflects the fact that at the collision energies of interest the nona- 
diabatic transition probabilities also vary by similar factors, and 
that many of them have extremely small values. These small 
transition probabilities are sensitive to variations in both the in- 
put quantum chemical data (potentials and couplings) and the 
method employed for solving the coupled channel equations in 
the quantum scattering calculationsQln addition, transitions can 
occur in symmetries other than for example, in the singlet 
system the p-, d-, and /-states can populate the ionic state via 
rotational couplings between n and £ states, as well as excite 
other states due to both rotational and radial nonadiabatic cou- 
plings in n, A, and other symmetries. The same is true for the 
triplet system with larger (at least by a factor of 3) initial statis- 
tical probabilities. For example, the cross-sections for 3s — > 3p 
excitation and their sensitivity to the input quantu m chemical 
data are discussed at length in Belvae v et alJ (1201 Oh . where we 
have shown that the results may vary by a factor of about 100. 
Similarly, from our investigations it is possible to make estimates 
of how much rate coefficients might vary for all transitions, and 
thus we define a "fluctuation factor" representing an estimate of 
this variation and which is somewhat analogous to an estimate 
of the uncertainty. The probable range for a given rate coefficient 
can be estimated to be from the value given in Table[T]multiplied 
by the minimum value of the fluctuation factor, up to the same 
value multiplied by the maximum value of the fluctuation fac- 
tor, noting that from the arguments discussed above, particularly 
contributions from other symmetries, it is generally more likely 
that the rate coefficient would be increased rather than decreased. 

We have assembled estimates of these fluctuation factors for 
all transitions in Table|2] A common feature of the quantum scat- 
tering calculations, particularly for very adiabatic transitions, is 
that the sensitivity of the cross-sections to the input quantum 
chemistry data usually decreases with increasing collision en- 
ergy, and thus cross-sections are often most uncertain very close 
to threshold, becoming less uncertain with increasing collision 
energy. This means that rate coefficients are usually most un- 
certain at low temperatures, becoming less uncertain at higher 
temperatures. Thus, in some cases we provide a range of fluctu- 
ation factors corresponding to the lower and higher ends of the 
temperature range considered. The rate coefficient for 3 s — > 3p 
excitation is roughly 10~ 21 cm 3 /s at 2000 K, somewhere in the 
middle of the range for all transitions. Thus, it is quite natural to 
expect that fluctuation factors for transitions with small transi- 
tion probabilties can be much larger than 100, while the fluctua- 
tion coefficients for the transitions determined by large transition 
probabilities should be close to unity. 

It is worth making clear that the smallest rate coefficients 
usually have the largest fluctuation factors: even if these rates 
were multiplied by such a large fluctuation factor they would 
not become dominant or substantially change the general pattern 
of relative strengths between different transitions. Nevertheless, 
for astrophysical applications it might be important to have such 
estimates of the possible uncertainties for the rate coefficients in 
order to see if any such transitions could have effects and thus 
warrant more accurate study. 



1 Note, in Belvae v" ;t al.l d2010h we found satisfactory agreement be- 
tween two different methods and codes when calculating cross-sections 
for the 3s — > 3p and As transitions. These transitions have reason- 
ably large cross-sections. However, recent additional calculations have 
shown some differences between calculations for transitions with much 
smaller cross-sections. 
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Table 1. Rate coefficients (crv), in units of cm 3 s _1 , for selected temperatures in the range T = 500-8000 K, for the excitation and 
deexcitation processes Na(nl)+H(ls)^>Na(n'l')+H(ls), or where indicated ion-pair production Na(n0+H(ls)— »Na + (2i 2 2p 6 )+H _ 
and mutual neutralisation Na + (2s 2 2/? 6 )+H~ ->Na(nT)+H(ls). 



Initial final state n'Y 



state nl 


3s 




3p 




4s 




3d 




4p 


5s 




4d 




4f 




5p 




Na + +ET 




















500 K 






















3s 






1.07E- 


-37 


1.09E- 


-50 


3.41E- 


-54 


8.22E-56 


2.75E- 


61 


8.05E- 


-64 


2.76E- 


-64 


9.36E- 


-66 


5 64E- 


■60 


- F 


5.81E- 


-17 






1.77E- 


-23 


5.45E- 


-28 


2.03E-29 


8.25E- 


■36 


9.86E- 


-40 


5.60E- 


-40 


3.41E- 


-43 


4 64E- 


■34 


4s 


1.65E- 


-18 


4.93E- 


12 






4.15E- 


-14 


7.63E-16 


2.77E- 


-22 


2.93E- 


-27 


2.61E- 


-27 


6.14E- 


-30 


4 14E- 


-20 


3d 


2.02E- 


-18 


5.95E- 


-13 


1.62E- 


-10 






6.19E-13 


3.91E- 


-20 


6.67E- 


-24 


2.21E- 


-24 


3.40E- 


-27 


2 20E- 


17 


4p 


1.92E- 


-18 


8.72E- 


-13 


1.17E- 


-10 


2.44E- 


-1 1 




1.42E- 


-14 


7.98E- 


-19 


4.21E- 


-19 


2.20E- 


-23 


2 53E- 


16 


5 s 


8.97E- 


-20 


4.96E- 


■15 


5.96E- 


-13 


2.15E- 


-14 


1.99E-10 






3.45E- 


-15 


8.97E- 


16 


2.01E- 


-18 


5 75E- 


-15 


4d 


2.60E- 


-21 


5.87E- 


-18 


6.24E- 


-17 


3.64E- 


-17 


1.10E-13 


3.42E- 


-14 






6.18E- 


-10 


1.83E- 


-16 


8.69E- 


■17 


4f 


7.34E- 


-22 


2.75E- 


18 


4.57E- 


-17 


9.91E- 


-18 


4.80E-14 


7.32E- 


-15 


5.09E- 


-10 






1.64E- 


-16 


5 80E- 


- 1 7 


5p 


2.08E- 


-22 


1.40E- 


-20 


9.02E- 


-19 


1.28E- 


-19 


2.10E-17 


1.37E- 


-16 


1.26E- 


-15 


1.38E- 


-15 






7 48E- 


■19 


Na + +H~ 


3.67E- 


-15 


5.55E- 


-10 


1.77E- 


-07 


2.41E- 


-08 


7.04E-09 
2000 K 


1.15E- 


11 


1.75E- 


-14 


1 .42E- 


-14 


2.18E- 


-17 






3 s 






1.71E- 


-21 


7.65E- 


26 


3.82E- 


-27 


7.85E-28 


7.45E- 


-30 


1.09E- 


-31 


8.46E- 


-32 


2.76E- 


-31 


5 51E- 


.27 


3d 


7.86E- 


-17 






4.35E- 


-15 


1.51E- 


-16 


4.25E-17 


2.63E- 


-20 


8.94E- 


-24 


8.19E- 


-24 


2.72E- 


-25 


5 87E- 


17 


4s 


8.48E- 


-18 


7.21E- 


-12 






7.36E- 


-1 1 


1.08E-11 


1.95E- 


-15 


1.93E- 


-19 


2.41E- 


-19 


4.63E- 


-21 


2 72E- 


1 j 


3d 


1.00E- 


-18 


5.91E- 


-13 


1.74E- 


-10 






7.55E-12 


1.50E- 


-15 


1.31E- 


-17 


8.22E- 


-18 


1.97E- 


-19 


8 30E- 


-12 


4p 


7.57E- 


-19 


6.12E- 


-13 


9.35E- 


11 


2.77E- 


1 1 




1.50E- 


11 


2.08E- 


-14 


2.31E- 


-14 


7.31E- 


-17 


9 04E- 


-12 


5 s 


1.78E- 


-19 


9.40E- 


-15 


4.19E- 


-13 


1.37E- 


-13 


3.73E-10 






3.71E- 


-13 


1.75E- 


-13 


1.67E- 


-15 


4 63E- 


.13 


4d 


1.39E- 


-21 


1.69E- 


-18 


2.20E- 


-17 


6.35E- 


-16 


2.74E-13 


1.97E- 


-13 






1.90E- 


-09 


3.43E- 


-13 


5.85E- 


16 


4f 


7.94E- 


-22 


1.15E- 


18 


2.04E- 


-17 


2.94E- 


-16 


2.25E-13 


6.86E- 


■14 


1.41E- 


-09 






2.08E- 


-13 


4 75E- 


16 


5p 


8.32E- 


-21 


1.22E- 


-19 


1.26E- 


-18 


2.26E- 


-17 


2.29E-15 


2.10E- 


■15 


8.16E- 


-13 


6.69E- 


-13 






7 33E- 


-18 


Na + +H" 


2.49E- 


-15 


3.96E- 


-10 


1.11E- 


-07 


1.43E- 


-08 


4.24E-09 
4000 K 


8.74E- 


-12 


2.08E- 


-14 


2.28E- 


-14 


1.10E- 


-16 






3s 






2.02E- 


■18 


2.28E- 


-21 


1.43E- 


-22 


3.76E-23 


7.97E- 


-25 


3.28E- 


-26 


3.09E- 


-26 


2.07E- 


-25 


2 61E- 


-21 


3p 

- F 


3.02E- 


-16 






1.51E- 


-13 


1.04E- 


-14 


3.50E-15 


5.73E- 


-18 


4.08E- 


-21 


3.32E- 


-21 


9.77E- 


-22 


4 48E- 


-14 


4s 


2.40E- 


-17 


1.07E- 


-11 






1.99E- 


-10 


3.79E-11 


1.49E- 


-14 


3.16E- 


-18 


3.29E- 


-18 


4.19E- 


-19 


7 03E- 


■10 


3d 


1.04E- 


-18 


5.03E- 


-13 


1.37E- 


-10 






9.22E-12 


4.43E- 


-14 


1.11E- 


-15 


6.69E- 


-16 


4.21E- 


-17 


6 05E- 


j j 


4p 


6.74E- 


-19 


4.20E- 


-13 


6.45E- 


11 


2.28E- 


-1 1 




7.51E- 


-1 1 


2.96E- 


-13 


4.15E- 


-13 


6.83E- 


-15 


4 49E- 


_l 1 


5 s 


1.23E- 


-19 


5.93E- 


-15 


2.19E- 


-13 


9.46E- 


-13 


6.48E-10 






4.01E- 


-12 


1.79E- 


-12 


3.54E- 


-14 


1 14E- 


-12 


4d 


1.65E- 


-21 


1.38E- 


-18 


1.51E- 


-17 


7.73E- 


-15 


8.31E-13 


1.31E- 


12 






2.81E- 


-09 


4.74E- 


-12 


1.53E- 


■15 


4f 


1.13E- 


-21 


8.12E- 


19 


1.14E- 


-17 


3.38E- 


-15 


8.48E-13 


4.24E- 


-13 


2.04E- 


-09 






3.00E- 


-12 


1 58E- 


1 s 


5p 


2.08E- 


-20 


6.56E- 


-19 


3.99E- 


-18 


5.83E- 


-16 


3.82E-14 


2.30E- 


-14 


9.44E- 


-12 


8.21E- 


-12 






6 78E- 


17 


Na + +H~ 


3.50E- 


-15 


4.03E- 


■10 


8.98E- 


-08 


1.12E- 


-08 


3.37E-09 
6000 K 


9.95E- 


12 


4.09E- 


-14 


5.79E- 


-14 


9.09E- 


-16 






3 s 






5.24E- 


■17 


1.25E- 


-19 


5.85E- 


-21 


1.52E-21 


4.13E- 


-23 


4.65E- 


-24 


4.76E- 


-24 


3.03E- 


-23 


3 22E- 


■19 


3p 

- F 


1.02E- 


-15 






5.55E- 


-13 


3.74E- 


-14 


1.27E-14 


2.77E- 


-17 


1.20E- 


-19 


1 .72E- 


-19 


4.01E- 


-20 


4 66E- 




4s 


6.02E- 


-17 


1.37E- 


-11 






2.37E- 


-10 


4.83E-11 


2.33E- 


-14 


2.33E- 


-17 


2.45E- 


-17 


5.12E- 


-18 


2 03E- 


()9 


3d 


1.28E- 


-18 


4.19E- 


-13 


1.08E- 


-10 






7.87E-12 


2.28E- 


-13 


6.94E- 


-15 


6.06E- 


-15 


5.61E- 


-16 


1 13E- 


■10 


4p 

^F 


7.24E- 


-19 


3.09E- 


-13 


4.78E- 


11 


1.71E- 


1 1 




1.56E- 


-10 


1.29E- 


-12 


1 .40E- 


-12 


5.44E- 


-14 


7 38E- 


. 1 1 


5 s 


1.19E- 


-19 


4.08E- 


-15 


1 .40E- 


-13 


3.00E- 


-12 


9.44E-10 






1.11E- 


-1 1 


5.94E- 


-12 


1.84E- 


-13 


1 69E- 


-12 


4d 


3.71E- 


-21 


4.91E- 


-18 


3.87E- 


-17 


2.53E- 


-14 


2.16E-12 


3.08E- 


12 






3.43E- 


-09 


1.68E- 


-11 


3.57E- 


-15 


4f 


2.75E- 


-21 


5.07E- 


■18 


2.94E- 


-17 


1.60E- 


-14 


1.70E-12 


1.19E- 


12 


2.48E- 


-09 






1.16E- 


-11 


3.10E- 


■15 


5p 


4.54E- 


-20 


3.07E- 


18 


1.59E- 


-17 


3.84E- 


-15 


1.71E-13 


9.56E- 


■14 


3.14E- 


-11 


3.01E- 


-11 






2.91E- 


16 


Na++IT 


6.22E- 


-15 


4.62E- 


10 


8.16E- 


-08 


9.94E- 


-09 


3.01E-09 
8000 K 


1.13E- 


-11 


8.66E- 


-14 


1.04E- 


-13 


3.76E- 


-15 






3s 






3.77E- 


16 


1.35E- 


-18 


5.00E- 


-20 


1.25E-20 


4.13E- 


-22 


1.01E- 


-22 


1.03E- 


-22 


5.00E- 


-22 


4.95E- 


-18 


3p 


2.66E- 


-15 






1.15E- 


-12 


6.84E- 


-14 


2.24E-14 


6.04E- 


■17 


1.13E- 


-18 


2.54E- 


-18 


6.87E- 


-19 


1.66E- 


-12 


4s 


1.39E- 


-16 


1.68E- 


11 






2.40E- 


-10 


4.77E-11 


2.80E- 


-14 


1.17E- 


-16 


1.79E- 


16 


4.34E- 


-17 


3.43E- 


■09 


3d 


1.90E- 


-18 


3.69E- 


-13 


8.88E- 


11 






7.07E-12 


6.16E- 


■13 


1.84E- 


-14 


2.49E- 


-14 


3.13E- 


-15 


1.51E- 


10 


4p 


9.70E- 


-19 


2.46E- 


■13 


3.60E- 


11 


1.44E- 


-11 




2.41E- 


-10 


3.05E- 


-12 


2.92E- 


-12 


1.91E- 


-13 


9.36E- 


-11 


5 s 


1.62E- 


-19 


3.36E- 


-15 


1.07E- 


-13 


6.36E- 


-12 


1.23E-09 






2.02E- 


-11 


1.28E- 


-11 


6.63E- 


-13 


2.13E- 


-12 


4d 


1.02E- 


-20 


1.61E- 


-17 


1.14E- 


-16 


4.86E- 


-14 


3.96E-12 


5.15E- 


-12 






3.91E- 


-09 


3.91E- 


-11 


6.39E- 


-15 


4f 


7.43E- 


-21 


2.60E- 


-17 


1.26E- 


-16 


4.74E- 


-14 


2.73E-12 


2.36E- 


-12 


2.82E- 


-09 






2.72E- 


-11 


5.04E- 


-15 


5p 
Na + +H~ 


9.13E- 


-20 


1.78E- 


-17 


7.73E- 


-17 


1.50E- 


-14 


4.51E-13 


3.08E- 


-13 


7.12E- 


-11 


6.88E- 


11 






8.48E- 


16 


1.15E- 


-14 


5.46E- 


10 


7.75E- 


08 


9.22E- 


-09 


2.81E-09 


1.26E- 


•11 


1.48E- 


-13 


1.62E- 


-13 


1.08E- 


-14 







These calculations have been employed to determine the to- 
tal rate coefficient for mutual neutralization. This rate coefficient 
can be approximated as 2.44 x 10~ 7 (300K/r) 032 cm 3 /s with an 
error of less than 5% for 500 <T< 80 00 K. This new result ca n 
be up to 50% larger than that given by Dickinso n et al.1 (jl999): 
the difference is primarily because a very much finer energy grid 
has been employed in the lBelvaev et alj d2010t) work. 



3. Discussion 

The rate coefficients show a general patter n similar to that found 
for Li+H collisions (Barklem et aD l2003h . with large rates for 
charge exchange processes, particularly from the first excited 
s state, and comparatively small rates for excitation processes, 
except in the cases of transitions between neighbouring excited 
levels. One may therefore reasonably expect similar effects on 
the statistical equilibrium calculations from inclusion of these 
processses: i.e. excitation and deexcitation processes have little 
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Table 2. Estimated maximum values of the "fluctuation factors" (a measure of the uncertainty) in the rate coefficients, given as 
factors of the rate coefficients in TableQ] The minimum value is estimated to be 0.5 in all cases. If a range is given, the first number 
corresponds to 2000 K and the second the 8000 K temperature result. For example, if the table lists 10-2, then the maximum 
fluctuation factor in the rate coefficient is estimated to be 10 near 2000 K and 2 near 8000 K. 



Initial final state n'l' 

state nl 3p 4s 3d 4p 5s 4d 4f 5p Na + +H~ 



3s 


120-70 50-10 10 


20 


100-10 


100-10 


100-10 


10 J -100 


2 


3p 


2 5 


2 


5-3 


100-10 


10 3 -30 


10 5 -100 


2 


4s 


6 


2 


100-10 


10 4 ~400 


10 4 -100 


10 5 -10 3 


2 


3d 




2 


20-2 


100-20 


10-2 


10 4 -2 


2-3 


4p 






2 


2 


2 


2 


2 


5s 








2 


2 


100-2 


10 


4d 










2 


10 3 -2 


10 4 -100 


4f 












10-2 


100^10 


5p 














10 4 -30 



10° 




or \ 

1 0" 25 l< " 

1Q -25 1Q -20 1Q -15 1Q -10 1Q -5 1Q 

Rate coefficient (this work) [cm 3 /s] 



Fig. 1. Comparison of rate coefficients at 6000 K from the 
Drawin formula (y-axis) with those from this work (x-axis). The 
dotted line is the one-to-one relation. The horizontal lines show 
the range in our values implied by the fluctuation factors, which 
are taken from Tableland log-interpolated at 6000 K between 
the extremes if a range is given. 



or negligible effect, while charge exchange processes have sig- 
nificant effects. It is, however, important that this be confirmed 
via detailed statistical equilibrium calculations. 

Regarding comparison of th ese data with those from the 
commonly used Drawin formula (|Steenbock & Holw eger 1984; 
lDrawinll968HDrawin & Emardll973l) . the situation is again sim- 
ilar to that for Li. The Drawin formula gives rate coefficients 
for optically allowed transitions (the Drawin formula cannot be 
used for optically forbidden transitions or charge exchange) be- 
tween one and seven orders of magnitude larger than the results 
presented here, see Fig. Q] Note that there is some correlation 
between the Drawin results and ours which simply stems from 
the fact that the Drawin formula correctly accounts for the fact 



that transitions with smaller energy thresholds will have larger 
rate coefficients. In any case there is a scatter of some six orders 
of magnitude around this correlation and a significant offset. We 
also note that the differences cannot be explained within the es- 
timated fluctuation factors for the results, and in addition there is 
no correlation between the difference and the uncertainty in our 
rate coefficient. 
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